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The recent discovery by the Daya-Bay and RENO experiments, that 613 is nonzero and relatively 
large, significantly impacts existing experiments and the planning of future facilities. In many 
scenarios, the nonzero value of #13 implies that #23 is likely to be different from n/4. Additionally, 
large detectors will be sensitive to matter effects on the oscillations of atmospheric neutrinos, making 
it possible to determine the neutrino mass hierarchy and the octant of 823- We show that a 50 kT 
magnetized liquid argon neutrino detector can ascertain the mass hierarchy with a significance larger 
than 4<r with moderate exposure times, and the octant at the level of 2 — 3<r with greater exposure. 
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Introduction: Neutrino oscillations in the standard 
three-flavor framework are parametrized by (a) two mass- 
squared differences Am^ = m 2 — to 2 , j — 2, 3; (b) three 
mixing angles 9ij\ and (c) one CP phase Sep- Over 
the past year, the T2K % MINOS @], and Double- 
Chooz S\ experiments have provided evidence that the 
mixing angle #13 is nonzero and not significantly smaller 
than the upper bound set by the Chooz [4] experiment. 
Recently, the Daya-Bay and RENO experiments have 
provided clinching evidence that #13 ^= at more than 
5cr: sin 2 29 13 = 0.089 ±0mi(stat)±0.005(syst) 5], and 
sin 2 20i3 = 0.113 ± 0.013(stai) ± 0.019(syst) [6], respec- 
tively. 

These results open new windows of opportunity for 
present-day or near-future neutrino experiments. If #13 
were tiny, an elaborate and expensive program would 
have been needed to measure the undetermined neutrino 
parameters. Its moderately large value allows for mea- 
surable matter effects [7[ in both beam and atmospheric 
neutrino experiments. This will yield better than an- 
ticipated precision in measurements of neutrino param- 
eters, and impacts the planning of future neutrino fa- 
cilities. Specifically, it allows for a credible case to be 
made for building advanced technology, large-mass at- 
mospheric neutrino detectors. The resolution of impor- 
tant outstanding questions, such as the mass hierarchy 
(i.e., whether TO3 > toi or 777,3 < r7 ^i) an d the octant of 
#23 (i.e., whether $23 is larger or smaller than 7r/4), may 
come within the purview of such detectors. 

The mass hierarchy plays a crucial role in the formula- 
tion of any theoretical effort designed to carry us beyond 
the standard model (SM) because models with a normal 
hierarchy (NH) are significantly and qualitatively differ- 
ent from those with an inverted hierarchy (IH); see, e.g., 
|3|. Thus, knowledge of mass hierarchy can help in dis- 
criminating between classes of models and consequently 



sharpen the focus of our search for new physics. A rela- 
tively large value of $13 suggests that #23 is unlikely to be 
maximal if the breaking of a /u — r exchange symmetry [9( 
in the lepton sector causes #13 to be nonzero. A nonzero 
value of 0i3 also opens the door to measurements of CP 
violation in the neutrino sector. 

Atmospheric neutrino detectors have been very impor- 
tant as discovery tools in the past. While they may not 
match the precision to pin down the energy and direc- 
tion of an event characteristic of long baseline experi- 
ments, they have some strong compensating advantages. 
Atmospheric neutrinos offer a broad range in baselines 
(L ~ 20 km to 12500 km) and energy E (100 McV 
to 10 TeV) that can be tapped into by such an exper- 
iment. An important consequence of this is the res- 
olution or alleviation of degeneracies endemic to long 
baseline beam experiments [10J. Large-mass liquid ar- 
gon time projection chambers (LATPCs) are based on 
one of the most promising technologies for charged par- 
ticle detection. They have unprecedented capabilities for 
the detection of neutrino interactions, rare events, and 
dark matter because of their precise and sensitive spatial 
and calorimetric resolution. The ICARUS [ll| and Ar- 
goNeuT [12] detectors have provided excellent examples 
of the capabilities of liquid argon as a neutrino interac- 
tion detector medium. 

In this Letter we study how well the hierarchy and oc- 
tant may be determined by a liquid argon detector using 
atmospheric neutrinos in light of the recently reported 
nonzero value of $13. 

Experimental specifications: We consider a large 



liquid argon detector as discussed in [13H15J . which can 
detect charged particles with good resolution over the 
energy range of MeV to multi-GeV. Magnetization over a 
50-100 kT volume has been proposed [16| and we assume 
it in what follows. We also assume the following detector 



resolutions over the GeV energy ranges relevant to our 
calculations |13l.[l7j: 

o\E e = ct_e m = 0.01 , 

°E had = v / (0-15) 2 /^ad + (0.03) 2 , 

C6> c = 0.03 radians = 1.72°, 

CT e M = V6 h «d = °- 04 radians = 2.29° , 

where E na d is the hadron energy in GeV, ge are the 
energy resolution widths in GeV and ag are the angular 
resolution widths of electrons, muons, and hadrons. 

Since E v = Ei ep + E na d, the neutrino energy resolution 
width, for both f„ and z/ e , is 



°E V 



V(0.01) 2 + (0.15) 2 /(y^) + (0.03) 2 , (1) 



where we have used Ehad = yE Vl y being the rapidity. 
In our computation, we take the average rapidity in the 
GeV energy region to be 0.45 for neutrinos and 0.3 for 
antineutrinos [18| . The angular resolution of the detec- 
tor for neutrinos can be worked out to be ag uc = 2.8°, 
ag v = 3.2° [13]. Charged lepton detection and separa- 
tion (e vs /x) without charge identification is possible for 
Eiepton > few MeV. The charge identification capability 
of the detector is incorporated as discussed in [l9[ . For 
electron events, we have conservatively assumed a 20% 
probability of charge identification in the energy range 
1 — 5 GeV, and none for events with energies higher than 
5 GeV. The muon charge identification capability of a 
LATPC is excellent for energies 1-10 GeV and we have 
assumed it to be 100% 

The atmospheric fluxes are taken from the three- 
dimensional calculation in [20] • The Earth matter profile 
defined in [21| is used to take into account matter effects 
on the oscillation probabilities. 

Hierarchy sensitivity: The hierarchy sensitivity 
from atmospheric neutrinos is computed with marginal- 
ization over the following test parameter ranges: 

• 6*23 from 38° - 52° (sin 2 23 = 0.38 - 0.62) , 

• lAmfi] from (2.05 - 2.75) x 10~ 3 eV 2 , 

• 013 from 5.5° - 11.0° (sin 2 20 i3 = 0.04 - 0.14) , 

• 5cp from — 2tt . 

We take the true values to be (lAm 2 ,^)^ = 2.4 x 10~ 3 
eV 2 , (Scp)tr = 0, (sin 2 23 ) tr . = 0.4,0.5, and 0.6 which 
allows for both maximal and nonmaximal values, and 
(013 )tr over the 3a range of recent measurements. The 
solar parameters Am|i and 0i2 are fixed at the values 
of 8 x 10~ 5 eV 2 and 34°, respectively, since the effect of 
their variation within uncertainties is negligible. The true 
value of Sep is chosen to be zero throughout these calcu- 
lations, since the principal contribution to both the hier- 
archy and the octant sensitivity is from the muon survival 
probability, which has a weak dependence on (5cp)tr- So 
the sensitivity is largely unaffected by the value of this 



parameter [22j . We assume the true hierarchy to be nor- 
mal and compute the ability of the detector to rule out 
the inverted hierarchy and vice versa. Flux uncertainties 
and systematic errors are included, using the method of 
pulls as described in [22j: flux normalization error 20%, 
flux tilt factor, zenith angle dependence uncertainty 5%, 
overall cross-section uncertainty 10%, and overall system- 
atic uncertainty 5%. The number of bins is chosen to be 
9 energy bins in the range 1 — 10 GeV and 18 cos0 z bins 
in the range —1.0 to —0.1. 

We also take into account the uncertainties of neu- 
trino parameters in the form of appropriate priors with 
the following la ranges: cr(|Am 31 |) = O.OSlAm 2 ^, 
cr(sin 2 20i 3 ) = 0.01, and cr(sin 2 20 23 ) = 0.02 sin 2 20 23 . 

The hierarchy sensitivity with priors \tot is obtained 

from a combined minimization of X 2 + X 2 + Xprior- 

Figure [1] depicts the hierarchy sensitivity for both NH 
and IH as a function of (0i 3 )t r for (sin 2 23 )t r = 0.4, 0.5, 
and 0.6. The Daya-Bay and RENO best-fit points are 
indicated in the figures. In general, the hierarchy sensi- 
tivity from both muon and electron events improves with 
an increase in the true value of 0i 3 and 23 . However, 
the sensitivity arising from muon events is considerably 
greater due to the statistical advantage enjoyed by muon 
events and the superior charge identification capability 
for muons as compared to electrons (22|, [23[ . For both hi- 
erarchy and octant calculation below, sensitivity strongly 
depends on matter effects, the detection of which relies 
on charge identification. Additionally, as also discussed 
in [l9|, 112, [23| , the excellent angular resolution of LATPC 
enables a full exploitation of the variations in the survival 
probability P w (the main constituent of muon events) 
with baseline. 

In Table I we list the \ 2 values for hierarchy sensitivity 
for (sin 2 20 13 ) tr = 0.1 and (sin 2 23 ) tr = 0.5. To highlight 
the role of the measurement uncertainty in 0i 3 , we pro- 
vide x 2 values for three cases: (a) no prior, (b) prior 
with cr(sin 2 20i 3 ) = 0.01 (present uncertainty), and (c) 
cr(sin 2 20i 3 ) = 0.005 (expected uncertainty in the near 
future). As expected, prior knowledge of sin 20 i3 has a 
dramatic effect on the ability to determine the hierarchy. 
However, an improvement in cr(sin 20i 3 ) below the cur- 
rent value does not give a concomitant improvement in 
the sensitivity. We note that without priors, the hierar- 
chy sensitivity is much greater for a true NH, since in this 
case the sensitivity arises from resonant matter effects in 
the muon events, while for a true IH the matter resonance 
occurs in the antimuon events. The atmospheric v^ flux 
is about twice the v^ flux, leading to greater sensitivity 
for NH. The prior terms tend to move the test parame- 
ters at the Xmin closer to their true values. This causes a 
tension between xft + xl and Xprior resulting in an overall 
increase in the total x 2 - If the true hierarchy is inverted, 
then the test parameter values at the Xmin without priors 
are further removed from the true values as compared to 
NH. Also, the sensitivity for sin 2 23 = 0.5 becomes su- 
perior to that for sin 2 23 = 0.6, as can be seen from the 
right panel of Fig. [TJ we have checked that without the 
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FIG. 1: Marginalized hierarchy sensitivity including priors in a liquid argon detector with 250 kT yr exposure as a function 
of (sin 2 29n)tr for (sin 2 $23)tr — 0.4,0.5, and 0.6. The vertical lines represent the Daya-Bay and RENO best-fit values. The 
left panel is for a true normal hierarchy (NH) and the right panel is for a true inverted hierarchy (IH). Sensitivities for the 
nonmagnetized version of the detector for sin 2 #23 = 0.5 are also shown. 
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25.6 


11.7 


0.01 


35.8 


35.4 


0.005 


36.1 


37.2 



TABLE I: Values of marginalized hierarchy sensitivity 

(x 2 tot P "° r ) f° r various values of er(sin 2S13) with 250 kT yr ex- 
posure for (sin 2 20i3)t r = 0.1 and (sin 2 6*23)tr = 0.5. The second 
column is for a true normal hierarchy and the third column is for 
a true inverted hierarchy. 



prior term the x 2 s for IH for different values of 6*23 follow 
the same behavior as that for NH. Clearly, the effect of 
priors is more significant for IH. 

For comparison we also present the sensitivity for 
a nonmagnetized version of the detector in Fig. [T] for 
sin 2 6*23 = 0.5. Note the crucial role played by the charge 
identification capability of a magnetized detector in dis- 
cerning the mass hierarchy. 

In summary, it is evident that if the value of #13 lies 
within the la range preferred by Daya-Bay /RENO, then 
with a 50 kT detector running for 5 years, an ~ 6c hi- 
erarchy determination is possible for both hierarchies for 
sin 2 26»i 3 = 0.1 and (sin 2 6»2 3 )tr = 0.5 (c/. Table I). 
Hence, even a reduced exposure of 100 kT yr may give 
a 4<t hierarchy discrimination irrespective of whether the 
hierarchy is normal or inverted. 

Octant sensitivity: For this analysis, we adopt the 
viewpoint that once the hierarchy is determined with 250 
kT yr exposure, the octant may be tackled with addi- 
tional exposure and a priori knowledge of the hierarchy. 
We assume a total exposure of 500 kT yr and consider 
the octant discrimination separately for NH and IH. As 
with hierarchy discrimination, we include priors and per- 
form a combined minimization of x 2 + x 2 + Xprior- The 



principal contribution to the octant sensitivity is again 
from muon events due to the sin #23 term in the survival 
probability (24J ]. 

In Fig. [31 we plot the Xmin values indicating the ability 
of the experiment to rule out the wrong octant for two 
values of #13 close to the Daya-Bay and RENO best fits. 
In the left panel of Fig. [2] the hierarchy is assumed to 
be normal and for the right panel it is inverted. 

We see that a 2a discrimination is possible for |6*23 — 
7r/4| > 3.5° for values of #13 close to the present best-fit, 
if the hierarchy is normal. That is, the octant discrimi- 
nation is possible only if the value of sin 2 26*23 is less than 
0.985. For larger values of #23, octant discrimination is 
difficult. For smaller values of sin 2 26*23, the value of Xmin 
increases rapidly and a 3<r octant discrimination is pos- 
sible for I6I23 - 7r/4| > 5° or sin 2 26* 2 3 < 0.97 for a normal 
hierarchy. If the hierarchy is inverted, octant sensitiv- 
ity is worse, and only a 2a discrimination is possible for 
1 023 - tt/4| > 4° or sin 2 2<9 23 < 0.98. 

Conclusions: The nature of the neutrino mass hierar- 
chy and the octant of 6*23 are vital to our efforts to build 
theories beyond the SM. With regard to these two issues, 
we have explored the implications of the recent Daya-Bay 
and RENO results on 6*13 for the planning of future exper- 
iments. In particular, we have demonstrated the excep- 
tional capability of a large-mass magnetized LATPC to 
determine the hierarchy to high significance with mod- 
erate exposure times. The detector is also sensitive to 
the octant of 6*23 although with a lower significance. Our 
results highlight the superior capability of a magnetized 
detector as compared to one without magnetization. 
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FIG. 2: Marginalized octant sensitivity including priors in a liquid argon detector with 500 kT yr exposure as a function of true 
#23 for (sin 2 29i3)tr = 0.07 and 0.1. The left panel is for NH and the right panel is for IH. 
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